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The goal of the paper is to determine boundary conditions in PDE models of collisions 
I of microswimmers in a viscous fluid. We consider two self-propelled spheres (microswim- 

• • mers) moving towards each other in viscous fluid. We first show that under commonly 

^ used no-slip boundary conditions on the fiuid-solid interface the microswimmers do not 

collide which is a generalization of the well-known no-collision paradox for solid bodies 
(with no self-propulsion) in a viscous fiuid. Secondly, we show that the microswimmers 
^ do collide when the no-slip boundary conditions are replaced by the Navier boundary 

conditions which therefore provides an adequate model of microswimmers such as swim- 
ming bacteria. The self-propulsion mechanism generates a drag force pulling a bacterium 
backwards and the collision problem is reduced to the analysis of competition between 
the drag and self-propulsion. For no-slip this is done by utilizing the Lorentz Reciprocal 
Theorem and the analytical solution for two solid spheres in the fluid. The analysis for the 
Navier boundary conditions is based on the variational formulation of Stokes problem. 
A Poincare type inequality for symmetrized gradient is introduced in this work. 

Keywords: Self-propulsion, Stokes flow, no collision paradox, Navier boundary conditions 

AMS Subject Classiflcation: 76D07,76D08 



1. Introduction 

In recent years active suspensions became a focus of attention in both biophysical 
and mathematical communities. This was motivated by a number of experiments 
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highlighting differences between active and passive suspensions and suggesting novel 
engineering applications. These experiments (just mention a few) include: decrease 
of effective viscosity for active suspensions compared with passive ones with the 
same volume fraction of inclusions ^ , swimming bacteria rotating small asymmetric 
gears ^, enhanced mixing in active suspensions of swimming bacteria ^, collective 
behavior of swimming microorganisms E. coli observed through correlations in 
swimmer orientations on scales order of magnitude larger than the size of individual 
swimmer. We refer also to theoretical works Swimmer-swimmer interactions 
are essential for all of the above phenomena. Therefore, for the mathematical models 
to recover the physically observed phenomena it is important to model accurately 
these interactions. 

Our focus will be on one of the aspects of swimmer-swimmer interactions - colli- 
sions. It is hard to perform an experiment with collisions of two swimmers. Instead, 
collisions between swimmers and walls have been clearly observed in physical ex- 
periments 8.9.10^ -^^g sJiow in this paper, some of the commonly accepted models 
of active swimmers in the fluid considered in literature do not allow for collisions 
in finite time. In this paper we propose a model that captures finite time collisions. 

Modeling of active suspensions consists of several components: fiuid motion, 
particle (swimmer) motion, particle-fluid interactions, and self-propulsion. Modeling 
of fiuid motion at the scale of microswimmers (micron scale) by Stokes equation 
has been universally accepted. Particles are commonly taken to be rigid bodies. The 
interactions between the particle and a fluid are prescribed through the boundary 
conditions on the surface of the particle. Existing models of self-propelled swimmers 
use no-slip boundary conditions on the entire swimmer's surface or its part''''^^"'^^''^'^, 
which mean that fluid sticks to the surface and moves with the same velocity as 
the boundary. Here we model self-propulsion by a point source term with given 
magnitude and direction applied at some distance behind the body'''^^. We present 
the model in details in Section 2. 

For passive suspensions (no self-propulsion) there is a well-known no-collision 
paradox It states that two rigid particles with no-slip boundary conditions, 
immersed in a Stokesian fluid and pushed towards each other with a constant force 
will take infinite time to collide. This is a consequence of an asymptotic relationship 
V ^ h between the velocity v of the spheres and the distance 2h between the spheres 
in the limit of small h, which leads to asymptotically exponential decay of the 
distance between the particles. It is called a paradox due to a mismatch between 
the theoretical predictions of no collisions and the experimental observations of 
possibility of collisions. 

A possible remedy to the no-collision paradox is to replace the no-slip boundary 
conditions by Navier ones. The Navier boundary conditions allow for a slip on the 
fiuid/solid interface with friction linearly proportional to the slip velocity. The pro- 
portionality coefficient ^ can be related to the surface roughness. Navier boundary 
conditions were derived as a rigorous homogenization limit of models of immersed 
rigid bodies with mean roughness (size of bumps) of order j3 when P is small (see 
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^^). Two passive spheres with Navicr conditions on the boundary pushed towards 
each other by a constant force will collide in finite time. This follows from an es- 
timate on the drag force for a two spheres moving towards each other with a unit 
velocity, which was recently proven in 

The no collision paradox also holds for active suspensions in the following sense. 
Two swimmers moving towards each other in Stokes fluid do not collide in finite 
time, see Section 4. Wc also prove that for sufficiently large initial velocities two 
swimmers with Navier boundary conditions, moving towards each other, collide in 
finite time, see Section 4. The no-slip boundary conditions have been successfully 
used in the study of dilute suspensions of nonintcracting swimmers"^®''''"'^®, pairwise 
interacting swimmers^^'^**'^" and collective behavior of swimmers^^'^. 

To recover collisions for more concentrated suspensions (where interactions are 
important) the no-slip conditions should be replaced by Navier boundary conditions 
(or, possibly by some other conditions). 

The key ingredients of our analysis are estimates for the drag force generated 
by the propulsion mechanism onto the body of the swimmer. In the analysis of the 
model with no-slip boundary conditions we use the Lorentz reciprocal theorem (or 
the second Green's formula for the Stokes equation, see Theorem 5.1) to reduce the 
original problem with singular ^-functions to an auxiliary problem without singu- 
larities. For the reduced problem the exact solution is known in a form of an infinite 
scries. The analysis of the model with the Navier boundary conditions is based on 
the variational principle for nonhomogeneous Stokes problem, the representation 
of the drag force via minimization principle (see Proposition 6.1), and the special 
Poincare-type inequality (Proposition 6.2) arc used. 

The paper is organized as follows. The models for swimmers with no-slip and 
Navier boundary conditions are presented in Section 2. In Section 3 we outline 
some relevant results for passive spheres. The main results (no collisions for no- 
slip boundary conditions, finite time collisions for Navier boundary conditions) are 
presented in Section 4. The proofs of these results are presented in Sections 5 and 
6, respectively. 

2. PDE/ODE model for microswimmers 

We consider the motion of two self-propelled spheres (swimmers). At time t>0 they 
are separated by distance 2h{t) and occupy regions Bi — |x : |x — xj| < l} and 
^2 = {x : |x - x2| < l}, see Fig I.Here xj = (0, 0, -l-h{t)) and x^ = (0, 0, l+h{t)) 
are centers of swimmers' bodies. Assume that the swimmers are immersed in a 
viscous incompressible fluid. The fluid occupies the domain 0,h{t) = ^^\{Bi U B2). 
For a given distance h > the fluid is described by a vector field u which solves the 
Stokes problem in Qf^: 

-Au + Vp = ^ 5(x - xi)fpd\ V • w = 0. (2.1) 

i=l,2 
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Here (5(x) is a delta- function, fp is a positive parameter which represents the inten- 
sity of self-propulsion, = (0, 0, 1) and = (0, 0,-1) are directions of the swim- 
mers' motion, = (0, 0, —1 — /i — A) and = (0, 0, 1 + /i + A) are ends of bacteria 

flagella and A is the length of flagella. We also assume that V{u — w$) e [L^(f]?i)] ^, 
where u$(x) = g{x-xl)d^ +g{x-xl)<P and ^(x) = l/87r(l/|x| • I + xx/|x|3) is 
the Oseen tensor. 

Equations (2.1) are supplemented by one of two following boundary conditions 
on dBii 

u + h' {t)d' = 0; (2.2) 
{u + h\t)d') ■ n = 0, {u + h\t)d')xn=-2l3[D{u)n]xn. (2.3) 

The boundary conditions (2.2) and (2.3) arc no-slip and Navicr boundary condi- 
tions, respectively. Here (3 is a positive parameter, and D{u) = ^(Vm -|- (Vu)'^). 
Note that if /3 = 0, then (2.3) becomes (2.2). The vector u + h'{t)d^ is the relative 
velocity of fluid particles near a swimmer with respect to this swimmer. Navier 
boundary conditions (2.3) state that tangential components of relative velocity and 
the stress on the boundary are proportional. 

For given flow u and pressure p the drag force on sphere dBi is 

^drag(w) = - j cr{u,p)nds ■ d\ (2.4) 

dBi 

where a{u,p) = —2D{u) + pi is the stress tensor and n is a normal vector on dBi. 
Due to the symmetry of the problem (spheres are identical) the right hand side of 
(2.4) does not depend on i = 1,2. The force balance equation is similar for both 
spheres and looks as follows 

mh"{t)-J^dr^^{u) + fp = 0, (2.5) 

where m > is mass of a swimmer. 

We assume that the swimmers do not rotate and, thus, torque balance equation 
trivially holds. 

Remark 2.1. In the balance equation (2.5) the propulsion enters explicitly in the 
third term, fp, and implicitly in the second term, J^irag(■w). We wish to rewrite 
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(2.5) in such a way that the propulsion enters only in one of the terms and does 
so explicitly. We can consider the drag force J\jrag(w) as a function of distance h{t) 
and speed h'{t) in the following sense. For a given distance h{t) and velocity h'{t) 
the function u is uniquely determined by the Stokes problem (2.1) with boundary 
conditions (2.2) or (2.3). Next our goal will be to explicitly write out the dependence 
of ^drag(w) on h{t) and h'{t). Due to the linearity of the Stokes problem we can 
decompose the fluid flow u into two components 

u = -h'{t)v + fpW, (2.6) 

where the first component —h'{t)v is generated by the motion of the spheres and 
the second component fpW is due to the propulsion force. The function v and w 
solve the following Stokes problems 



' Av - Vpi) = in Qm), 

V • i; = in n^^t^ , 

(w - d') • n = on dBi, i = l,2, 
S{v-d') xn = -2/3 [D{v)n] x n, 



Aw - Vp^ = 5{x- xP)d' in n^^f^, 

i=l,2 

V • w = in 

w • n = on dBi, i = 1,2, 
w xn = -2/3 [D{w)n] x n. 



Due to the linear dependence of the drag force Jdrag on the flow u we obtain 

^ drag('^^) h (t)Kpags ^ fp^'propj 

where Kpass '■= -7\irag(i') is the drag force on a passive sphere moving with a unit 
velocity, and Kprop := .7^drag(w) is the drag force on a motionless sphere due to a 
unit propulsion force. Note that both coefficients Kprop and Kpass depend on the 
distance h. 

Using the drag coefficients Kprop and Kpass) we can rewrite the balance equation 
(2.5) as follows: 

mh"{t) + Kp,,,h'(t) + fp{l - Kp,op) = 0. (2.7) 

Now, the self-propulsion parameter enters only the third term in (2.7), unlike 
(2.5) where J\irag(M) depends on fp. If fp = 0, then (2.7) becomes a force balance 
equation for passive spheres. 



3. Review of known results about collision for peissive inclusions 

We now briefly review the result for passive spheres. 

Consider two unit spheres moving along towards each other along the common 
axis z (see Fig l.).Let h{t) be a half-distance between the spheres. Then the speed 
of each sphere is —h'{t). Assume that the external force /oxt pushes the spheres 
towards each other. Denote by Tdiag a magnitude of the drag force of the fluid 
computed by the formula (2.4). The force balance on each sphere is 



^ drag ~t~ /ext — 0* 



(3.1) 
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Prom remark 2.1 the drag force J^drag depends hnearly on the speed of the 
spheres —h'{t): 

drag — l^passh (t). ("^•2) 

The drag coefficient Kpass is a function of half-distance h. It also depends on the 
parameter /? > 0, i.e., it depends on the choice of the boundary conditions. 

The dependence of Kpass on h for no-sHp boimdary conditions, /3 = 0, is known 
as a classical result (see, e.g., ^^'^^). For Navier boundary conditions, /3 > 0, the 
result has been justified recently (see ^'^): 

1, 1 
1 

h' 

First, consider no-slip boundary conditions. From (3.4) with /3 = we may 
conclude that for some positive constant C > the following inequality holds: 

C 

/^pass < ~J^- (>^-5) 

Using (3.1), (3.2), and (3.5), we obtain 

hit) > /i(0)exp(-C/extt). 

It imphes that h{t) can not vanish for any finite t > 0. It says that the distance 
between spheres will never become zero, thus they will never collide whatever pos- 
itive numbers h{0) and /ext are. This contradicts to what happens in real life and 
it is known Uke no-collision paradox. 

Secondly, consider Navier boundary conditions. It is sufficient to restrict our- 
selves to the case h < (3 « 1. The relation (3.3) implies that for some positive 
constant C > the following inequality holds: 

C 1 

Kpass > 'p^^^fl' (^-6) 

Using (3.1), (3.2) and (3.6), we obtain 

Ch{t){\nh{t) - 1) - Ch{0){\nh(0) - 1) > /exti- 

This inequality together with h'{t) < predicts that there exists Tcoii such that 
< Tcoii < 00 and h{Tcoii) = 0. Thus, in the case of Navier boundary conditions 
collisions do occur. 

4. Two main results 

4.1. No collisions for swimmers with no-slip boundary conditions 

The following theorem is the main result for the problem with the no-slip boundary 
conditions. 
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Theorem 4.1. Consider coupled PDE/ODE model (2.1), (2.2) and (2.5). Let h{t) 
be a half-distance between two swimmers at t > 0, introduced in Section 3. Then 
there exist positive constants Ci > and C2 > such that 



h{t) > Cie-^=* 



(4.1) 



for all t>0. 



I.: ( 



' Av - Vpy 
V-v = 0, 
V = Vd\ 

K drag('^) 



0, 



Proof. For simplicity, consider the case without inertia (m = 0). The case m > 
relies on the same arguments. Let u satisfy (2.1), (2. 2) and (2.5). Then u = v + w 
where v and w satisfy the following systems of equations: 

' Aw-Vp^= E ^(x-xf)/pd% 

i=l,2 

IL: I V • w = 0, 

w = Wd\ 

, •^drag(w) = 0. 

In both systems, I and II, equalities in the first and the second line are satisfied 
in , equalities in the third line are satisfied on dBi for all i = 1,2, equalities in 
the fourth line are satisfied for all i — 1,2. li fp and h are given we can find v and 
q from the first, the second and the fourth equalities of system I. After that we 
can find ^ G M from the third equality. Similarly, for system II. Thus, there exist 
mappings Gy and defined by problem I and problem II, respectively, such that 
Qv{h, fp) = V and Qw{h, fp) = W (like the Stokes law for a sphere J\irag = QniiU). 
The result follows from the following two propositions: 

Proposition 4.1. Let V G M. be the number obtained from system I (i.e., V = 
Qv{h, fp) ). Then there exists a positive constant C > such that 

V < Cfph. (4.2) 

The result of proposition 4.1 follows from the relations (3.5) (see ^^-^i) and 

fp ~ ^'^pass* 

Proposition 4.2. Let W G M. be the number obtained in system II (i.e., W = 
Qw{h, fp) ). Then the following inequality holds: 



W <0. 



(4.3) 



The proof of proposition 4.2 is relagated to Section 5. 

Now (4.2), (4.3) a.ndu = v + w imply that -h'{t) = V + W < Cfph{t). Substi- 
tuting it to (2.5) we have 



h'{t) > -Cfph{t), 



which implies the statement of the theorem. □ 
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Fig. 2. Dependence of ftprop on flagella length A with h = 0.01 and /3 = 0.1 

4.2. Finite time collisions for swimmers with Navier boundary 
conditions 

Theorem 4.2. Consider coupled PDE/ODE model (2.1), (2.3), (2.5) and m > 0. 
Then there exists a ■positive constant C > which may depend on h{0), m, (3 and 
A such that if h'{0) > C, then there exits time Tcoii > such that h{Tcoii) — 0. 

Tcoii is time of collision. Qualitative etimates on Tcoii which can be compared 
with measurements have been obtained in (see also Remark 4.2). The proof 
of Theorem 4.2 is relagated to Subsection 6.2. The proof is based on variational 
principle and the Poincare type inequality (see Subsection 6.1). 

Remark 4.1. The Theorem 4.2 states that swimmers collide if the initial speed is 
sufficiently large. Thus, the result on passive spheres is extended on active swim- 
mers. But one can expect that swimmers must collide for arbitrary initial speed and 
separation distance. This is due to the propulsion force pushing swimmers towards 
each other which is absent in the case of passive spheres. In order to prove collisions 
without any restrictions on initial conditions and m > one need to show 

1-Kprop > ;(A,/3), (4.4) 

where Kprop was introduced in Section 3 and represents the drag force on a swimmer 
generated by self-propulsion = 1 of both swimmers. The function l{X,f3) is a 
strictly positive function, it may depend on A and /?, but does not depend on the 
swimmer's separation h. The inequality (4.4) means that the drag force due to 
propulsion mechanism ((5-functions in the Stokes equation) is not greater than the 
actual propulsion. In other words, the pushing mechanism really pushes swimmers 
towards each other. 

For fixed roughness parameter /3 = 0.1 and small separation distance h = 0.01 
the force Kprop is a function of tail length A, denote it by g{X)- Numerics suggest 
that inequality (4.4) holds for all A, /3 > 0. It is optimal in the following sense: when 
A the function g{X) = Kprop 1 (see Fig 2). This makes the inequality (4.4) 
difficult to prove. 

Remark 4.2. It is also interesting to estimate the time before collision Tcoii in 
terms of parameter /3. Time before collision Tcoii can be evaluated by the following 
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formula 

rHO) 

T'coll 



where U{h) = —h'{t). The following estimate for the collision time holds (/ext = 
1, h{Q) = 1): 

Teoii~l//31n(l//i). (4.6) 

To estibhsh the formula (4.5), one can use (2.7), the assumption (4.4) and the 
relation for (3.3). Numerics on collision time are performed in our paper 



5. Proof of proposition 4.2 

Wc use the following well-known theorem 

Theorem 5.1 (Lorentz Reciprocal Theorem). Assume that equalities Au^ — 
Vp^ = -Fi, Av? - Vp^ = -F2 and V • = 0, V • = hold in domain Q.. Then 

j Fi ■ u^dx - j a{u^,p^)n ■ u^ds = j ' '^^^^ " j o■(w^3^^)n • u^ds. (5.1) 
a an n an 

Here n is an inward (for fl ) normal vector. 

Let us take Q, = ^u, = w, p^ = pw and -F2 = — ^ (5(x — xp/pd*. Let 

i=l,2 

be such function that Am^ — Vp^ = (in other words, Fi =0), and satisfies the 
following boundary conditions: 

= Wd' on dBi, i = 1,2. 

Prom (5.1) we get: 

i=l,2 g'g i=l,2 

The following equality holds: 

j 0'(u\p^)n = J"drag(w^)d', 

dBi 

and due to symmetry of the problem for [u^ ,p^) we have 

• u\yil) = d^ . u\xl) = -uli^l), i = 1,2, 
where u^(Xp) is a ^-component of u^(Xp). Thus, we get 

w = -j-\ry.{4). (5.2) 

dragV"' } 

Thus, it remains to prove positivity of Uz{^l). 
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Lemma 5.1. The following inequality holds 

uli^D > 0. (5.3) 

Proof. We recall that v} satisfies the following system 

f Au^ - Vjji = 0, in Vlh 
\ V • = 0, mQ.h 

v} = WA\ on dBi, i = 1,2. 

The proof is divided into three parts. First, the result on p. 247 is is used to 
write the exact solution of system above. Second, we calculate at a point x^, 
and, finally, a necessary estimate is obtained. 

STEP 1. Exact solution. Introduce cylindrical coordinates for x = {x,y,z) 

X = pcoscj), 
y = psincl), 
z = z. 

Denote = (cos sin 0) and = (0,0,1). Due to the axial symmetry of the 
problem we have = Upep + UzCz where scalar functions Up and Uz do not depend 
on (f). Following we introduce the stream function t/j = \p[p, z) 

Uz{z,p) = -^^V'(2,P), Up{z,p) = ^^V'(2,P)- 

In order to write function ip we introduce the bipolar coordinates C € [0, +oo), 
p, £ [0, tt] (we need the case z > only) 

sinh C sin r? 

z = c — — , p = c- 



cosh ^ — cos 7] cosh ( — cos r] 

where c = sinh a and a are such positive number that cosh a = 1 + h. Note that 
due to the following equality 

{z — ccoth^)^ + = (c cosech^)^ 

one can easily verify that surface {( = a} is the sphere dBi. Also we will need the 
formula: 

C + ^^ = lnf "Tt^'! V (5-4) 



^p + i{z- c) 

The function tjj may be written as follows (see ^4,23^ 

oo 

i;{C,7i) = (coshC - cos r,)-3/2^t/„(C)C-^f (COST?), 



where 



n=0 



1 3 
UniO = K sinh(n - -)C + rf„ cosh(n + -)C, 



^-1/2. X _ Pn-l{x) - Pn+l{x) 
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Here Pn{x) are the Legendrc polynomials 

Pn{x) = 



(.x^ - 1)" 



2"n! dx" 

and the coefficients 6„ and d„ are given by the following formulas 



K = 



WW'smh.^a^ — x 



"v^(2n- 1) 

4cosh^(n +\)a + 2{2n + 1) sinh^ a 
2 sinh(2n + l)a - (2n + 1) sinh 2a 



- 1 



dn = 



V2{2n + 3) 
4cosh^(n + - 2(2n+ 1) sinh^ a 
~ 2 sinh(2n + l)a - (2n + 1) sinh 2a 



STEP 2. The calculation of u^{x.p). Wc need to calculate when zq = 1 + h + ^ 
and p — >■ 0+. In bipolar coordinates it corresponds to 

Co = In ^° ^ and rio = 0. 
zo-c 

Then using (5.4) (to obtain ^ and P„(l) = 1, P^{1) = n(n + 1)/2 and the 
equality 

dtp dC dt/j di] dtp 
pdp pdp dC, dp pdrj 

we obtain 

1 d 



Uz{zo,0) 



P^o+ p dp 



,(cosh-iCo-l)-'/'E^"(Co) 



(5.5) 



n=0 



Since — h{2 + h) we have that — Zg < and the statement of proposition 
follows if we prove that ?7„(Co) > for all n e N. 

STEP 3. Positivity o/U'„(Co). The nonpenetration condition (u — d') ■n = implies 
the following relation between coefficients 6„ and d„ (see, e.g., formula (42) in ^^): 



sinh(m + l)a 
smh(m — l)a 



(5.6) 



where to = n + 1/2 and 

^ _ c2 to2-1/4 



(to + l)e-(™-i)" - (to - l)e-('"+i)" 



sinh(TO — 1)q! 



2\/2 to2 - 1 
The formula for Un{^) has the following form: 

sinh(TO + 1)^ sinh(TO + l)Q: 



sinh(TO — 1)^ sinh(TO — l)a 



sinh(TO - 1)^. (5.7) 
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Note that 

(m + l)e-('"-^)" - (m - l)e-(™+i)" = 2e"""(msinha + cosha) > 

Thus, 

Gm > 0. (5.8) 

Introduce /(^) := sinh(TO + 1)^ sinh(m — l)a — sinh(T7i — 1)^ sinh(TO + l)a. Consider 

^ G [0, a]. One can see that /(O) = /(a) = 0. Consider now M := max /(^) = 

£e[o,a] 

/(^o)- Assume that M > 0. Then 

f'i^o) = (to + 1)^ sinh(m + 1)^ sinh(TO — l)a — (to — 1)^ sinh(TO — 1)^ sinh(m + l)a 

> (m - 1)V(6) > 0, 

what is impossiple for maximum. It means that /(^) < for ^ e (0, a). Thus, 

sinh(TO + 1)^ sinh(TO + l)a ^ 
sinh(m — 1)^ sinh(m — l)a 

Combining (5.8) and (5.9) we obtain from (5.7) that J7„(^o) > 0. 

6. Proof of the collision result for Navier BC 
6.1. Aiixiliary statements 

Consider the following Stokes problem for smooth force density F: 



' -Au + Vp = F, 

V • M = 0, 

(m - d*) • n = 0, 
[ (w - d*) X n = -2/3 [D{u)n] x n. 



(6.1) 



Define energy functional: 

£^{v) = 2 J \D{v)\' + /^^ 1^ - - / (6-2) 

defined on the following class of functions: 

A={v\{v-d')-n = on dB'} . 
One can easily prove the following statement: 

Proposition 6.1. Let u solve Stokes problem (6.1). Then the following statements 
hold true: 

(i) u is a minimizant of on class A: 

min£'F(u) = £f{u). 
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C^'i-- 










• 







Fig. 3. Domains ft* and ft. 

(a) Drag force J-{u) — J^^i a{u,p)n ■ d' may be evaluated by the following 
formula: 

F{u) =£'f(u) + j Fu. 
We will also need the following Poincare type inequality: 
Proposition 6.2. Let v be the function such that t; • n = on dBi. Then 

\v\^<C j \D{v)\+ j \v\\ (6.3) 

where fJ* and U, are bounded domains such that f2* C 17 C ^h, and the set T :— 
dfl n dBi is not empty (see Fig. 3). 

Proof. 

For the sake of simplicity we consider the case when £7 and $7* are symmetric with 
respect to plane xy: if x = (a;, y, x) e W (or f2), then x*'^™ :— {x, y, —z) e W (or (l). 
STEP 1. Assume by contradiction that there exists such sequence 



that 



a. 



{Vn ■ Vn ■ TO. ^ Q On P} 

|D(v„)| ^0 and j |w„ 



0. 



(6.4) 



Let = { 7 + A X x| 7, A e M"^} be the linear subspace of L^(r2) ("rigid motion 

r - 1 3 
space") and H : L'^{Vl) 



— )■ y is the corresponding projection operator. Introduce 

(6.5) 



Wn := (/ — n)-y„. Then for all n there exit 7, A e E'^ such that 

w„(x) = 7„ + A„ X X + w„(x). 
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STEP 2. We claim that 



Indeed, the functions t«„ satisiy the Korn inequality (see, e.g., Theorem 2.5 in ^^): 

[\Vwn?<C [\D{w^)\\ (6.7) 



L2(fi) 



to V and, thus, they 



On the other hand, functions arc orthogonal in 
are orthogonal to constant three-dimensional vectors. It means that the avarage of 
each component of Wn in il is zero and Poincaret inequality implies that 

K|2<C / |V(«;„)r (6.8) 
Ja 

Collecting togeather (6.7), (6.8) and D{wn) = D{vn) wc get (6.6). 
STEP 3. We claim that there exists a sequence {n^ c N}^f^ 

as /c ^ +00. (6.9) 

Indeed, F is a finite dimensional space equipped with i^-norm. The sequence 

{7™ + An X X = t;„ - Wn] 



-ynk and 



is bounded in 



L2(0) 



n3 



and, thus, there exists a convergent subsequence 
+ X X ^- 7 + A X X. 



From /p |^'n|^ ^ and that T contains peace of the sphere dBi we get 



L 



I7 + A X xp = => 7 + A X X = for all X e r 



(6.10) 



and since span of all vectors x G F coincides with we get that 7 and A are zero 
vectors. This contradicts to the relation: 



= 1. 



(6.11) 

□ 



Remark 6.1. The inequality (6.3) without the integral over F 

j \v\'' <C J \D{v)\ (6.12) 

is not valid for function v such that w • n = on F. One can take nonzero vector 
A e {0} and the sequence: 

w„ = nA X (x — xl) 

Then w„ • n = u„ • (x - xji) = nA x (x - xj) • (x - xj) = 0, J^, |u„p = J^^ |A x 
(x-Xc)p -J> +00, while \D{vn)\ = 0. 
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Corollary 6.1. For functions v: [v — d') • n = the following holds 
j \v\^ <C j \D{v)\ +C j\v- df + C< C£{v) + C. 



(6.13) 



6.2. Proof of Theorem 4.2 

We consider the following Stokes problem 

-A« + Vp = Ei=i,2 /pdM(x - x;), 
V-u = 0, 

(ii + /i'(i)d*)-n = 0, 
[u + /i'(i)d') X n = -213 [D{u)-a\ x n, 

and the balance equation 

m/l"(t)- J-drag(M)+./p=0. 

where J-drag(it) = d* • 2D{u)n — pn. The definition of J-drag does not depend 
on i = 1, 2. 

STEP 1. First, we reduce the original system with singular force density to the 
problem with regular one. 

Denote K := —l/fph'{t). Consider the function u$ defined as follows: 



(6.14) 



[XWi\er H ^ [xmlez- 



where x ' 



r dz r dr 

[0, 1] is C°°-function such that 



X(x) 



-{I 



o,xeWA/2(x;). 



and t/j^ir, z) iRr ^ 
K 



V'4> 



1 1 

^{z - IpY + r2 ~ ^{z + IpY + r^_ 

where Zp = /i + 2 + Aisa coordinate of a tail. Function solves the following 
problem: 



i=l,2 

V • ?i = 0. 



(6.15) 



u$ • n = 0, 

u$ X n = -2/3 [£)(u$)n] x n. 

for some scalar and C°° function F which support belongs to Q* = 
{^/3 < |x — Xi| < ^/2,i = 1,2}. Consider u = —{1/K)u — {($. It satisfies the fol- 
lowing problem 

f-Au + Vp = F(x), 
V • w = 0, 
(m- d*) • n = 0, 
, (u - d*) X n = -2/3 [D{u)n\ x n. 



(6.16) 



December 27, 2012 2:52 WSPC/INSTRUCTION FILE 



Collisions-m3as'3 



16 L. Berlyand at el. 

STEP 2. We will prove the following inequality: 

J"drag(u) <C(|ln/l|+if2). 



Assume for simplicity K > 1. 
1. We claim that 



where 



/ 



Fu 



< -£{u) + C. 
4 



Indeed, in view of (6.13) 

/- 



Fu 



< KC 



(j^ j < KCV£{u) + 1 < ^£{u) + CK^ 



2. Consider £f{v) defined as follows 



Then from 1. we have that 



£f{u) > -£{u) - CK"^ 



Observe that 



^drag(w) = £f{u) + j Fu< £f{u) + ^£{u) + CK'^ 

< C£f{u) + CK^ < C min £{v) + CK'^ 



< C\\nh\+CK^. 

And thus inequality (6.17) is proven. 
STEP 3. The balance equation is 

mh"it) - J'drag(u) + fp = 0. 

The definition of J^rag does not depend on i = 1,2. 
Due to Step 2 the following inequality holds: 



mh"{t) - Ch'{t) In hit) + C/h'{t) + fp < 



(6.17) 



(6.18) 



(6.19) 



(6.20) 



£Fiv) = 2 J \D{v)\'' + IJ^- -Jf^ (6-21) 



(6.22) 

(6.23) 

(6.24) 
(6.25) 

(6.26) 

(6.27) 



Consider any v > 0. Assume that h'{t) < —u for all t e (0,t*) and /i'(t*) = —v. 
Then 



mh"{t) - Ch'{t) \nh{t) < -fp - C/h'{t) < -fp + Cjv 



(6.28) 
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Then by the integration 



mh'{t) < i-fp + C/u)t + mh'{Q) 

+C {h{t){\nh{t) - h{t)) - ;i(0)(ln/i(0) - /i(0))} 
< (-/p + C/i/)t + m/i'(0). 



Let us take —mh' 



(0) > C/vT + m{v + h{0) /T) then 



mh'{t) < -fpt - C/v{T -t)-mv- mh{0)/T. 



(6.29) 



Thus t^>T and -h'{0) < -h{0)/T. Thus, time of coUision T^oU is less then T. 
This completes the proof of theorem 4.2. 
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